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| Abstract 

We discuss quite surprising properties of the one-parameter family of modular 
(Auberson and Sabatier (1994)) nonlinear Schrodinger equations. We develop a 
unified theoretical framework for this family. Special attention is paid to the 
emergent dual time evolution scenarios which, albeit running in the real time 
parameter of the pertinent nonlinear equation, in each considered case, may be 
mapped among each other by means of a suitable analytic continuation in time 
procedure. This dynamical duality is characteristic for non-dissipative quantum 
motions and their dissipative (diffusion- type processes) partners, and naturally 
extends to classical motions in confining and scattering potentials. 
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An inspiration for the present paper comes from the recent publication p] discussing 

effects of various scale transformations upon the free Schrodinger picture dynamics. In 
particular, it has been noticed that an appropriate definition of the scale covariance 
induces Hamiltonians which mix, with a pertinent scale exponent as a hyperbolic ro- 
tation angle, an original free quantum dynamics with its free dissipative counterpart 
(effectively, a suitable version of the free Brownian motion), pQ c.f. also [2]. 

The two disparately different time evolution patterns do run with respect to the 
same real time label. However, the ultimate "mixing" effect of the above mentioned 
scale transformations takes the form of the the Lorentz-like transformation Eq. (35) 
where one Hamiltonian takes the role of a regular "time" while another of the "space" 
dimension labels. This obvious affinity with the Euclidean "space-time" notion and the 
involved complex analysis methods (e. g. Wick rotation, imaginary time transforma- 
tion, analytic continuation in time) sets both conceptual and possibly phenomenological 
obstacles/prospects pertaining to an existence of such dynamical patterns in nature. 
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Most puzzling for us is the apparently dual notion of the time label which while a 
priori referring to the real time evolution, may as well be interpreted as a Euclidean 
(imaginary time) evolution. It is our aim to address an issue in its full generality, by 
resorting to a one-parameter family of modular Schrodinger equations, where external 
conservative potentials admitted. 

The adopted perspective relies on standard approaches to nonlinear dynamical sys- 
tems where a sensitive dependence on a control (here, coupling strength) parameter 
may arise, possibly inducing global changes of properties of solutions to the equations 
of motion. That is exactly the case in the present analysis. We wish to demonstrate 
that a duality property (realized by an imaginary time transformation) does relate 
solutions of the modular nonlinear Schrodinger equation for various coupling constant 
regimes. 

In Section 2 we set a general Lagrangian and Hamiltonian framework for the sub- 
sequent discussion and indicate that effectively the dynamics can be reduced to three 
specific coupling value choices, [31 0] , each of them being separately discussed in the lit- 
erature. We aim at a stationary action principle formulation of the modular Schrodinger 
dynamics and the related hydrodynamics (the special cases are: the familiar Bohm- 
type quantum hydrodynamics [!5] and the hydrodynamical picture of diffusion-type 
processes [2]). Basic principles of the action principle workings are patterned after the 
classical hydrodynamics treatises, [HE]. 

In Section 3, in Hamiltonian dynamics terms, we give a new derivation and further 
generalize the original arguments of pQ. We employ scaling properties of the Shannon 
entropy of a continuous probability density p = if)*ip. The scale covariant patterns of 
evolution, non-dissipative with an admixture of a dissipative component, are estab- 
lished in external potential fields. 

Section 4 is devoted to a detailed analysis of the emergent time duality notion for 
the three special values 0,1,2 of the coupling parameter. There, we pay more atten- 
tion to a specific intertwine between confining and scattering dynamical systems. A 
properly addressed sign issue for the potential function appears to be vital for a math- 
ematical consistency of the formalism (links with the theory of dynamical semigroups 
that underlies the dissipative dynamics scenarios). As a byproduct of a discussion we 
establish the Lyapunov functionals for the considered dynamical patterns of behavior 
and the (Shannnon) entropy production time rate is singled out. 

Section 5 addresses more specific problems that allow to grasp the duality concept 
from a perspective of the theory of classical conservative systems and diffusion-type 
(specifically - Smoluchowski) stochastic processes. A number of illustrative examples 
is worked out in detail, with an emphasis on the time duality notion in the classically 
inspired (standard Hamilton- Jacobi equations) and dissipative patterns of evolution 
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(modified Hamilton- Jacobi: e.g. Smoluchowski diffusion processes vs quantum dynam- 
ics). 



2 Modular Schrodinger equations 



Let us consider a subclass of so-called modular Schrodinger equations [31 0] which are 
local and homogeneous nonlinear generalizations of the standard Schrodinger equation: 
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where ip(x,t) is a complex function, \ip\ = (ip*ip) 1 / 2 , V(x) is a real function and a 
coupling parameter k is non-negative, k > . 

If k > 0, the pertinent nonlinear dynamics is known to preserve the L 2 (R n ) norm of 
any initially given if), but not the Hilbert space scalar product (if), <f>) of two different, 
initially given if) and 0. The induced dynamics is non-unitary in L 2 (R n ); unitarity is 
restored if k = 0. The pertinent set of solutions is a subset of L 2 (R n ), but not a linear 
subspace (no superposition principle for k^O). 

For all k > a usual continuity equation holds true with the familiar quantum 
mechanical definition of a (probability) density current: d t p = —V ■ j, where p = if>*if) 
and j = (h/2mi)(if)*'Vif) — if)S/ip*). We consider normalized solutions only, which sets 
a standard form of j = p ■ v , where v = (h/2mi)[{Wif)/if)) — (Vif)* /if)*)] = (l/m)Vs is 
regarded as a gradient velocity field and p(x,t) = \if)\ 2 (x, t) is a probability density on 



2.1 Lagrangian formalism 

The one-parameter family of modular equations (pQ), together with its complex con- 
jugate version, derives from the local Lagrangian density £ by means of the sta- 
tionary action principle [H]-[5]. Let us consider a functional of ^-functions, their 
space and time derivatives, including complex conjugates: I[ip,ip*} = L(t)dt where 
L(t) = f C(x,t) dx, (we leave unspecified, possibly infinite, integration volume). 

We impose the stationary action condition 5I[if),if)*] = for independent variations 
8if), 5ip* which are bound to vanish at integration volume boundaries. By invoking 
elements of the pedestrian functional calculus [5], like e. g. 5C/5ip = dC/dif) — 

V 'i[dC/ 'd(y V0)], one ends up with the Euler-Lagrange equations: d t [dC/d(d t if)*)] = 
8C/8if)* and d t {d£/d(d t ip)\ = 8£/8ip. 

If we properly specify the Lagrangian density C = C K : 

c K (x, t) = y [v w) - - |^vv> • vr - v(x) if) r+ (2) 
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the stationary action principle yields a pair of adjoint modular equations which com- 
prise Eq. (pQ) in conjunction with its complex conjugate: 
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We have previously introduced a current velocity field v(x, t) by means of j = p ■ v. 
Our gradient assumption v = (l/m)Vs follows from an implicit reference to the familiar 
Madelung substitution: ip = \ip\ exp(is/h), where \ip\ 2 = p. Its consequent exploitation 
allows us to rewrite the Lagrangian density (j2J) as follows (we leave intact an original 
order of respective entries): 



C K (x, t) = -p d t s + y (u 2 + v 2 ) + V(x) -k^u 2 



(4) 



where u(x,t) = (h/2m) Vp/p is another velocity field (named an osmotic velocity). 
Here, 5I[p, s] — gives rise to a continuity equation d t p = — V(p • v) and yields: 



d t s + — (Vs) 2 + V+(1-k)Q = 
2m 



(5) 



where, in view of = p 1//2 , the familiar notion of the de Broglie-Bohm quantum 
potential, [HUH], naturally appears: 

2" 



h 2 Ap 1 ' 2 
Q = Q(x,t) = -— = -— 



4m 



Ap _ 1 /Vp 
~J~2\p~ 



The modular Schrodinger equation (1) takes the form: 

ihd t i> = \-(h 2 /2m)A + V}ip-KQip. 



(6) 



(7) 



2.2 Hamiltonian formalism 

A symplectic structure can be associated with the dynamical system (l)-(3) by in- 
troducing fields and n^* that are conjugate to ip and ip* respectively: = 
dC/d(d t i/j) = (ih/2)il)* and n^* = dC/d(d t i/j*) = —(ih/2)ijj. The subsequent Legendre- 
type transformation defines the Hamiltonian density: 



h 2 

H K = n^- dti> + 7i>* • d t ip* - C K = — W> ■ Vip* + 

2m 



P 



171 2 T T /-I \ 171 2 

— v 2 + V + (1 - k) — u 2 



h 2 fVtp* 

8m V tp* ip 



7i s d t s- £„ 



2 ' 2 

where, this time with respect to the polar fields p(x,t) and s(x,t), we have: n p 
dC/d(d t p) = and n s = dC/d(d t s) = —p. 
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The Hamiltonian reads H K (t) = J 7i K (x,t) dx. For the variational calculus it is not 
C but L — J Cdx that really matters. Therefore, it is useful to note that 

L(t) = -(d t s) - H K (t) , (9) 

where, in view of J pdx = 1, we can introduce the mean value (dts) = J pdtsdx. 

Let us evaluate the mean value of the generalized Hamilton- Jacobi equation (J5J). By 
assuming a proper behavior of p at the integration volume boundaries [TT], we readily 
get (Q) = fQpdx = +(m/2)(u 2 ) > 0. Thence, on dynamically admitted fields p{t) 
and s(x, t), L(t) = 0, i. e. (d t s) = —H K . 

Let us consider two function(al)s A = J A(x, t) dx and B = J B(x, t) dx, which may 
explicitly depend on time t. We define their Poisson bracket 

{A > B} = -Hj dx {wW-W*w-) (10) 

In particular, while identifying A = ip(x,t) and B = H K (t), we get the modular 
Schrodinger equation ([1]) in the form 

d t i, = ty, H K } (11) 

while, by setting A = i/)*, an adjoint equation arises 

a t r = {r, h k }. (12) 

We recall e. g. that 'k^ = —5H K /5ip while ip = 5H K /8^. 

Since the time dependence of H K (t) is realized only through the canonical fields, 
the Hamiltonian surely is a constant of motion. Thence (dts) as well. 

The polar decomposition ip = p 1 ^ 2 exp(is/H), ip* = p 1//2 exp(— is/K) preserves a 
symplectic structure. In the self-explanatory notation there holds 

{A, B} = {A, B}^ r = {A, B} p , s (13) 

and thence: 

r tt i 

d tP = {p, H K } = - 1 f = --V ( P Vs) (14) 

while 

d t s = {s, H K } = = --^(Vs) 2 - V - (1 - k) Q . (15) 

op 2m 

The result is valid for all k > 0. 

Let G(t) = f dxQ(x, t). Clearly, c.f. [6], if the time-dependence of G is realized only 
through canonically conjugate fields p(x,t) and s(x,t) (and their derivatives), then 

^§ = {G,H K }. (16) 

One should realize that a particular time-dependence pattern of G(t) critically relies 
on the chosen parameter range for k G R + . 
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2.3 Reduction to effective k, = 0, 1 and 2 self-coupling regimes 



In the notation of Eq. (TjJ, a distinguished role of k = 0, 1 or 2 coupling parameter values 
is particularly conspicuous. The relevance of only these three k — 0, 1, 2 parameter 
values may be further enhanced. 

Namely, let us recall [3] important differences between properties of solutions of 
Eqs. ([I]) and ([3]), depending on whether 0</c<l,«=lor/c>l. 

(i) In the parameter range < k < 1 (in fact < k < 1), if ip(x,t) = exp(is/h) 
actually is a solution of (1), then ip'(x', t') = \ip'\ exp(is'/h), with x' = x, t' = (1 — K) l ^ 2 t, 
\ip'\(x',t') = \ip\(x, (1 - kY^H') and s'(x',t') = (1 - k)- 1 / 2 s(x, t), automatically solves 
the linear Schrodinger equation: 

h 2 



ihdtiip' 



2m 



-A 



1 



V 



K 



(17) 



The scaling transformation replaces a nonlinear problem by the linear one, albeit with 
a re-scaled potential. 

(ii) For a specific value k = 1 we encounter the formalism that derives (though 
actually generalizing) from the wave picture of classical Newtonian mechanics, [5]. 

(iii) In case of k > 1, a repetition of previous scaling steps, provided we replace 
(1 — k) 1 / 2 by (k — l) 1 / 2 in the pertinent formulas, results in the following outcome: 
ip'(x',t') = exp(is'/h) is a solution of the nonlinear Schrodinger equation 



ihd[il)' 



& A 

A + - 

2m k 



V 



^' + 2 



h 2 A\ip'\ 
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Effectively, in the whole coupling parameter range k g R + , only the cases of k — 0, 
/t = 1 and k = 2 form a mutually exclusive family, both on mathematical and physical 
grounds. Since for k = and k = 2, the above scaling transformations trivialize, while 
being irrelevant for the distinctively "borderline" case of k — 1, we can safely restore 
the notation of Eqs. ((H) and (J3j). 

In connection with the choice of k — 2, one more observation is of utmost impor- 
tance, [3], see also [12]- [TBj. Namely, if a complex function 

(19) 



ip(x,t) = \ip\ exp(is / fr) 

is a solution of ([1]) with k = 2, then the reaZ function 

9*(x,t) = 1*01 exp(— s/fr) 

is a solution of the generalized (forward) heat equation 

ft 2 



2m 



A + V 



ft-! 



(20) 



(21) 
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with a diffusion coefficient D = h/2m and an external potential V = V/2mD. By 
setting V = we would arrive at the standard heat equation d t 9* = D A 6*. 

Another real function 0(x,t) = \ip\ exp(+s/h) is a solution of the time-adjoint 
(backwards) version of Eq. (12TT) : 



- hd t 6 



r ft 2 

2m 



e. 



(22) 



Note that if we would have started from Eqs. (|2ip and fl22l) with the purpose to 
arrive at the modular equation ([1]) with k = 2, the ill-posed Cauchy problem would 
possibly become a serious obstacle. That in view of the occurrence of the backwards 
parabolic equation. 

This ill-posedness might be healed by invoking the theory of strongly continuous 
dynamical semigroups, [T21 EH H7j. To this end we need to choose V(x) to be a 
continuous function that is bounded from above, so that V = —V becomes bounded 
from below. Then Eq. (12T]) would acquire a "canonical" form of the forward diffusion- 
type equation related to the contractive semigroup operator exp(—Ht/ft): 



me 



-H9 ± 



2m 



V 



0* 



(23) 



Eq. (J23J), together with its time adjoint 

ftd t 9 = H0 = 



h 2 

2m 







(24) 



stand for principal dynamical equations of so-called Euclidean quantum mechanics [13] . 
while falling into a broader framework of the Schrodinger boundary data and stochastic 
interpolation problem, [121 [13] see e -g- also [H]- 

Note, that if one assumes that the external potential V'(x) is a continuous func- 
tion that is bounded from below, then the operator H = —(ft 2 /2m) A + V is (essen- 
tially) self-adjoint in L 2 (R n ). So, we have consistently defined unitary transformations 
exp(—iHt/h) in L 2 (R n ) and iftd t ip = Hip, together with plus —iftd t ip* = Hip*, (k = 
case), as their local manifestations. 

A careful analysis reveals [T2l [T3] that the standard Schrodinger picture dynamics 
can be mapped into the contractive semigroup (generalized diffusion equation) dynam- 
ics (23) and (24), by means of an analytic continuation in time. This mapping we 
reproduce in the book-keeping (Wick rotation) form; (it) — > t. 

We emphasize that one should not mystify the emergent "imaginary time" trans- 
formation. The time label t pertains to fairly standard (unequivocally real) dynamics 
scenarios. However the detailed patterns of temporal behavior for two motion scenarios 
are very different. 



7 



Nonetheless, various (real) functionals of dynamical variables, evolving in accor- 
dance with the chosen motion rule (say, Schrodinger's ip(t) andtp*(t), may be consis- 
tently transformed into the affiliated functionals that evolve according to another, we 
call it dual, (dissipative 8*(t) and 0(t)) motion rule. The reverse transformation works 
as well. We shall address this issue in below, from varied perspectives. 

3 Scale covariant patterns of evolution 

Let [Al] stand for an arbitrary (albeit fixed for the present purpose) unit of length. For 
a continuous probability density p on R 1 we can introduce its (dimensionless) Shannon 
entropy functional, also named differential entropy: [2] 

S(p) = -Jd L^j ([Al]p) ln([Al]p) = - f dxp H[Al]p) . (25) 

If p(x, t) depends on time, then S(p) = S(t) may evolve in time as well. 

The time rate equation for S(t) is devoid of any [Al] input and for all k > has 
the very same functional form: 

DS = D{S,H K } = -(uv) (26) 

Obviously, (•) stands for the mean value with respect to a probability density p(x,t). 

The above time rate formula does not depend on a specific unit of length, that is 
present in the definition of the dimensionless Shannon entropy. However, the entropy 
functional itself is sensitive to scaling transformations. 

Setting x' = x//3, where (3 > is the scale parameter, we have 1 = j p(x) dx = 
J ft p(f3x') dx'. Accordingly, the scale transformation x — > x' = x/ (3 induces a transfor- 
mation of the probability density in question: 

p(x)^ftp(ftx') = p'(x'). (27) 

Given the Shannon entropy of the density p(x), we can compare the outcome with 
that for the density p'(x'): 

S'{p') = - J dx' p'{x')\n([Al]p'{x')) = - J dxp(x)\n(P [Al]p{x)) = S{p) -ln/3. 

(28) 

Consequently, the x — > x/f3 scaling is equivalent to the sole change of the length unit 
[Al] — > f3 [Al] in the Shannon entropy definition. 

Since (3 is a fixed scaling parameter, the time rate formula (126]) is scale independent. 
Note that by setting (3 = expa, we get S'(p') = S(p) — a. 
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We observe that u'(x') = (3u{x) and the scale invariance of DS = —(uv) tells us 
that v'(x') = (l/j3)v(x). Since we assume v(x,t) to be the gradient field, v = (l/m)Vs, 
we readily arrive at the corresponding scaling property of s(x,t): 

s'(x') = j 2 s{x) — > v'(x') = lv(x) . (29) 

We demand furthermore that the (k = 1) equations d t p = — V(p • v) and d t s + 
(l/2m)(Vs) 2 + V = are scale invariant, e.g. they retain their form after introducing 
x', p'(x',t) and s'(x',t) instead x, p(x,t) and s(x, t) respectively. This implies an 
induced scaling property of V(x) — > V'(x'): 

V'(x') = L V {x) (30) 

and thus H[ = (l/p 2 )^. 
In view of 

Q'{x',t) =p 2 Q(x,t) (31) 

the general evolution equations for k / 1 are not scale invariant. Let us consider a 
cumulative effect of the above scaling rules upon the Hamiltonian 

H K = j ' dxp[(mv 2 /2) + V + (1 - n){mu 2 /2)} (32) 

according to: 

H K = H K (t) - H'Jt) = J dx'H(p',s')(x',t) . (33) 

Following the guess of [lj, originally analyzed in the case of V = 0, we take j3 = 
exp(a/2); there follows 

dx p [— v 2 + V] + exp(+a) / dxp(\ — K)—u 2 (34) 
which can be recast as the hyperbolic transformation with a hyperbolic angle a 

H' K (a) = cosh a H K — sinh a K K . (35) 
Here, in addition to H Q , we encounter a new Hamiltonian generator K K 



K K = J dxp[^v 2 + V-{1-k)^u 2 \ , (36) 

The difference between H K and K K is encoded in the sole sign inversion of the last 
(1 — n)(mu 2 /2) entry. 

We can as well analyze an effect of the scaling transformation upon K K) [T]: 

K' R (a) = — sinh aH K + cosh aK K . (37) 
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Hyperbolic rotations that are explicit in equations (15S]) and (157)) form a conspicuous 
Lorentz-type transformation (hyperbolic rotation of coordinates in Minkowski space) 
of the direct analogue (H K , K K , 0, 0) of the familiar Minkowski space vectors (p , p, 0, 0) 
and/or (ct, x, 0, 0). 

We recall that the time label t is left untouched by hitherto considered scale trans- 
formations. Consequently, like H K , the generator K K and the induced a-family of 
generators H' K (a), K' n (a) are legitimate Hamiltonian generators of diverse time evo- 
lution scenarios, all running with respect to the same time variable t. The pertinent 
motions arise through common for all Hamiltonians, a priori prescribed, symplectic 
structure Eq. (1T5)) and the generic time evolution rule (1TB)) . 

Remark 1: By re-tracing back the passage from the modular equation to the 
corresponding Lagrangian and Hamiltonian, we realize that K K can be associated with 
another (k — 2)-family (k > 0) of modular Schrodinger equations: 

ihd T i) = [-(h 2 /2m)A + V]ip + (« - 2)Q V • ( 38 ) 

to be compared with the originally introduced K-family 

ihd t il)= [-(h 2 /2m)A + V]^-KQip. (39) 

We shall not be elaborate on the general k > 2 parameter regime and, in view of 
the previously established reduction procedure (Section 1.3), we confine our further 
discussion to specific k = 0, 1 or 2 cases. 

Remark 2: Let us notice that Kq = H2, H\ = K\ and K2 = Hq. We have 
H[ = exp(—a)H[ and the following hyperbolic transformation properties hold true for 
H Q and H 2 : 

H' = cosh aH — sinha^ (40) 

and 

H' 2 = — sinh aH Q + cosh aH 2 . (41) 

These transformation rules are a generalization of those presented for the V = case 
by L. Brenig, [TJ. 

Remark 3: Hamiltonians of the form 

H = (m/2)(v 2 + u 2 ) + (V) (42) 

and 

H 2 = (m/2)(v 2 -u 2 ) + (V) (43) 

are known in the literature, [HJ H3] , and are interpreted to set quantum-mechanical and 
dissipative-dynamical frameworks respectively. Some elementary hints in this connec- 
tion can also be found in (TJ [15], HSJ E] • Coming back to our observation that a rescaling 
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of the dimensional unit [A/] — > expa[A/] induces a transformation S'(p') = S(p) — a 
of the Shannon entropy S(p) for a continuous probability density p, we realize that the 
choice of a < implies a sharpening of the resolution unit, hence an effective growth of 
the Shannon entropy with the lowering of a — — \a\. Indeed, we have S'(p') = S(p) + \a\. 
In this situation, effective (rescaled form) Hamiltonians H' and H' 2 , always have an 
admixture of both non-dissipative and dissipative components, H and H 2 respectively 



4 Dynamical duality 



Remembering that for general functionals G(t) of p and s, we have G = {G,H K }, let 
us consider a product J-{x, t) = —p(x,t) s(x,t) of conjugate fields s and n s = —p. The 
time evolution of 

F(t) = / dxT{x,t) = ~{s) 



looks quite interesting: 



dx 



s(x, t) 



6H. 
Ss 



p{x,t) 



5p J 



(44) 



(45) 



— dx p 



m 



-V — V — (1 — K 



m 



u 



L 2 v ' 2 

A new Hamiltonian-type functional has emerged on the right-hand-side of the dynam- 
ical identity (1451) . Let us denote 



dx p 



m 



m 



—v z ±V±(l-K) — u' 



(46) 



We point out that in contrast to previous scaling-induced Hamiltonians, the negative 
sign has been generated both with respect to terms (m/2)(u 2 ) and (V). The +V — > — V 
mapping was painfully lacking in the previous discussion to have properly implemented, 
the implicit, analytic continuation in time procedure. 
The previous motion rule rewrites as 



dF 
~dT 



{F,H+} = -H^(t) 



(47) 



where H£ = H is the time evolution generator Eq. (jHJ). Note that H£ is here a constant 
of motion, while H~ (t) is not. 

After accounting for the Poisson bracket ffl3|) , we encounter a complementary rela- 
tionship for the time evolution that is generated by the induced Hamiltonian H~: 



dF 
~dt 



{F 1 H-} = -Ht(t). 
Presently, H£ is a constant of motion, while H~ (t) no longer is. 



(4* 
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For each value of k G R + , we thus arrive at dual time-evolution scenarios, generated 
by Hamiltonians and H~ , respectively. The duality notion stays in conformity with 
our previous discussion (Section 1.3) of the "imaginary time" transformation. 

These dual motions, even if started from the very same initial t — data and 
allowed to continue indefinitely, do result in disparately different patterns of behavior: 
non-dissipative and dissipative, respectively. 

Nonetheless, if we admit that the a priori chosen generator of motion is then 
we are interested as well in: 

j2p j tt r 

— = -{H-(t),Ht} = = +2 J pv V[V + (1 - k)Q) dx (49) 

while, presuming that H~ actually is the evolution generator, in the dual evolution 
formula 

j2p j tt+ r 

= ~{H+(t), H-} = = -2 J p v V[V + (1 - k)Q] dx . (50) 

The right-hand-sides of Eqs. (|49p . fl50l associate the mean power transfer rates 
(gain, loss or none) to the dual Hamiltonian evolutions of F(t), generated by H£ and 
H~ respectively. This remains divorced from the fact that H+(t) and H~{r) actually 
are constants of the pertinent dual motions. 

For clarity of our discussion we shall confine further attention to L 2 (R n ) solutions of 
(1) and (3). We assume that the external potential V(x) is a continuous function that 
is bounded from below. If the energy operator H = —(h 2 /2m)A + V is self-adjoint, 
then we have consistently defined unitary transformations exp(—iHt/h) in L 2 (R) so 
that ihdtip = Hip holds true (« = 0). 

The case of k = 1 we associate with two classes of external potentials ±V(x), with 
+V(x) bounded from below. This will allow us to discriminate between the confining 
and scattering regimes. The "borderline" meaning of k = 1 can be read out from 
Eq. (J3J), where mu 2 /2 contributions cancel away. 

In conformity with our previous discussion of generalized heat equations, related to 
solutions of Eq. (TjQ) with k — 2, given +V(x) and H, we pass to a pair of time-adjoint 
parabolic equations: hd t 9* = —H9* and hd t 9 = HQ. 

Here, 9*(x,t) = [exp(— Ht/K) #*](x, 0) represents a forward dynamical semigroup 
evolution, while 6(x,T — t) = exp(+Ht/h) 6(x,T) stands for a backward one. Both 
are unambiguously defined in a finite time interval [0,T], provided one has prescribed 
suitable end-point data [12] . 

The corresponding modular Schrodinger equations (plus their complex conjugate 
versions) read: 

(i) K = ihd t ip = [-(h 2 /2m)A + V)ip 
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(ii) K = i = > ifid t i) = [-(fr 2 /2m)A ± V - Q}^ 

(iii) K = 2 =>■ z£<9 t V> = [-(h 2 /2m)A -V - 2Q]ip. 

The associated Lagrangian densities (@J in the (p, s) -represent at ion and the induced 
dynamical rules are worth listing as well. In addition to the continuity equation d t p = 
— V(p • v) we have valid the Hamilton- Jacobi type equations: 

(i) K = 0; £ = -p[d t s + {m/2){v 2 + u 2 ) + V] <9 t s + (l/2m)(Vs) 2 + (V + Q) = 

(ii) k = 1; C = -p[d t s + (m/2)v 2 ± V] d t s + (l/2m)(Vs) 2 ± V = 

(iii) k = 2; £ = -p [d t s + (m/2)(v 2 - u 2 ) -V]=> d t s + (l/2m)(Vs) 2 -(V + Q) = 0. 
On dynamically admitted fields p(t) and s(x,t), L(t) = 0, i. e. (dts) = —H. The 

respective Hamiltonians (jSJ) do follow: 

(i) H + = J dxp[(m/2)v 2 + V + (m/2)w 2 ] 

{ii)Hf l =jdxp[{m/2)v 2 ±V} 

(iii) £T- = / dxp [(m/2)t» 2 - V - (m/2)w 2 ] 

We emphasize that, from the start, V(x) is chosen to be a continuous and bounded 
from below function. In the definition of the above Hamiltonians there is no k la- 
bel anymore and a subscript "cl" refers to the classically motivated (Hamilton- Jacobi 
theory) wave formalism. 

The evolution equations for F = — (s), c.f. jTJ and [TH [2], clearly define dual pairs: 



and 



= {F, H+} = - 
= {F, H-} = - 

F = {F H+} 
F = {F,H-}- 



dx p 
dx p 



m 



2 
in 
~2 



V 



m 



u 



.... 2 Tr m 2 
— v 2 + V-\ u 



-H-(t), 
-H+(t) 



dx p 



dx p 



~m 
— 1 
. 2 

m 



r 



t; 2 + y 



(51) 
(52) 

(53) 
(54) 



It is instructive to notice that the functional F(t) in Eqs. (1521) and (|54p may consis- 
tently play the role of a Lyapunov functional, indicating the preferred sense of time 
("time arrow") in the course of the evolution process. Namely, if we take (V) > 0, 
the right-hand-side expression is negative definite. Hence F(t) is a monotonically de- 
caying function of time which is a standard signature of a dissipation process, c.f. the 
Helmholtz free energy and the relative entropy discussion for diffusion-type processes, 

[S3 El H2|. 

In the notation of section 1, c.f. Eq. (jSJ), we have H ± = and = Hf. The 
corresponding (dual) time rate formulas f|49|) . fl50|) do follow. 

The motion rules for F(t) can be given more transparent form by reintroducing 
constants of the respective motions. Then 

F(t) = -m{v 2 )(t) + H ± (55) 
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and 

F(t) = -m(v 2 )(t) + H± . (56) 

Here, the non-negative term m(v 2 )(t) should receive due attention, because of its ut- 
most importance in the study of the Shannon entropy dynamics, [HI EJ [19], where 
it represents an entropy production time rate. The latter is generated solely by the 
dynamical processes which are intrinsic to the system and does not involve any en- 
ergy/heat flow, in or out of the potentially dissipative system (that would need the 
notion of an external to the system, thermal reservoir). 

Since H + and H~ are constants of respective motions, F(t)—t are monotonically 
decreasing in time quantities. This property extends to the generated dynamics as 
well. 

The speed (slowing down, or acceleration), with which the above decay process may 
occur, relies on the specific dynamical pattern of behavior of (v 2 ). That is quantified 
by —md(v 2 )/dt, hence: 

F(t) = {F{t), H±} = ±2 J pvV{V + Q)dx . (57) 

For H d generated motions, we have: 

F(t) = {F(t), H±} = ±2 J pvVVdx . (58) 

We point out that a major distinction between the dual dynamical rules is encoded 
in the right-hand-sides of the above equations: the conspicuous sign inversion is worth 
contemplation. The related integrals have a clear meaning of the power transfer (re- 
lease, absorption or possibly none), in the mean, that is induced by time evolution of 
the pertinent dynamical system, [16j [2] . 

Let us stress that the "imaginary time" transformation, even if not quite explicit 
in our discussion, hereby has been extended to dynamical models of purely classical 
provenance. The pertinent Wick rotation connects confined and scattering motions, 
admitted to occur in a continuous and bounded from below (confining) potential +V(x) 
and in its inverted (scattering) counterpart —V(x), respectively. An obvious example 
of an inverted harmonic oscillator [20J is worth mentioning at this point. Parabolic 
potential barriers [21] and repulsive 1/r 2 potentials [22] belong to the same category. 

In the above discussion, the sign inversion issue is manifested in second time deriva- 
tives of various functionals. In view of ±W presence, we can identify this behavior as 
a remnant of the standard classical (Newtonian) reasoning, [12] : if the sign looks wrong 
with respect to the classical Newton equation (e. g.we have +VV), we can correct this 
"defect" by interpreting time t as an "imaginary time" it (or ir to avoid a notational 
confusion) . 
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5 Physics-related implementations of the dual dy- 
namics patterns: a (dis)illusion of an "imaginary 
time" 

5.1 Generalities 

With the notational conventions D = h/2m, b(x,t) = v(x,t) + u(x,t), while imposing 
suitable boundary conditions (e.g. p, bp, vp vanishing at integration boundaries or 
infinities), we can write the Shannon entropy time rate of change in a number of 
equivalent ways, [2j. 

DS = D{S,H K } = -(uv) = (v 2 ) - (bv) (59) 

and 

DS = D(Vv) = (Vb) + D(u 2 } . (60) 

The non-negative entry (1/D)(v 2 ) is interpreted as the entropy production rate in the 
considered dynamical system [T6J, El [19] . 
We note that 

D(u 2 ) = -D{Vu) = -(Q) (61) 

m 

so that the mean divergence of an osmotic velocity is always negative. Here (Q) > 
holds true for all finite times, [2]. The value can be achieved, if at all, only in the 
asymptotic regime t — > oo. 

To give a flavor of the time duality (and specifically the "imaginary time transfor- 
mation") connection, let us mention that the free Brownian motion can be embedded 
in the above scheme by setting 6 = and regarding D as the diffusion constant. We 
have v = —u = —DVlnp and Shannon entropy time rate takes the form of the de 
Bruijn identity DS = (v 2 ) = (u 2 ), [2j. 

The Shannon entropy S(t) of the Brownian motion grows monotonically in time, 
solely due to the entropy production {1/ D)(v 2 ){t). The latter, in turn, is known to be 
a decreasing function of time (at least in the large time asymptotic) and ultimately is 
bound to vanish at t — > oo. 

In case of the diffusion process in a conservative potential it is not Shannon entropy, 
but the Helmholtz free energy that takes the role of the Lyapunov functional and sets 
the "time arrow", c.f. Eqs. (1551 and see [2j[T9]. Note that F(t) —tH~ is a monotonically 
decaying in time quantity. 

For comparison, the Schrodinger picture quantum dynamics typically involves b ^ 0. 
In the special case of the free motion S > 0, hence S(t) grows indefinitely, [16]. In 
the large time asymptotic S — > 0, while the entropy production (1/D)(v 2 ) remains 
untamed and never vanishes while approaching a finite positive value. 



15 



However, the would-be natural property S(t) > is not generic for the quantum 
motion in external potentials, [23] • Nonetheless F(t)—tH + does monotonically decrease 
with time t — > oo indicating the Lyapunov functional-induced "arrow of time". Even, 
though this dynamics is manifestly non-dissipative. 

5.2 Harmonic oscillator and its inverted partner 

Let us consider a standard classical harmonic oscillator problem, where 

n = L + ^ v (62) 

is an obvious constant of motion for the Newtonian system p = mq = —mu 2 q, 

<l(t) — Qo cos + — — sin ut (63) 
mu 

pit) = po cos ut — muqo sin ut . 

Clearly H = p^/2m + (mu; 2 /2)gg is a positive constant. 

Let us perform an analytic continuation in time, by considering the Wick rotation 
t — > —it, paralleled by the transformation of initial momentum data po — > —ipo- Once 
inserted to the above harmonic oscillator expressions, we arrive at 

H-ipo = + (rnu 2 /2)q 2 = -H (64) 



and 



together with 



Q-ipni—it) = q(t) = go cosh cut — ^-sinhut (65) 



P-i P0 {—it) = +ip(t) = —ipo cosh cut + imujq Q sinh ut , (66) 
which simply rewrites as 

pit) = — p cosh cut + muqo sinh ut . (67) 

We observe that 

q(—t) = q cosh out + sinh ut (68) 
mu! 

and 

— p(—t) = po cosh ut + muq sinh ut (69) 

are the familiar inverted oscillator solutions, generated by H, [2D] . 

Indeed, equations of motion for q(t) and p(t) directly derive from the Hamiltonian 
H-i PQ = —H with 

¥= £-^ 2 ' 2 (70) 
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They give rise to the Newton equation p = mq = +muo 2 q. However, the dynamics 
generated by H is related to that generated by —H by the time reflection: the latter 
dynamics runs backwards, if the former runs forward. 

Remark: At the first glance, the harmonic oscillator example may be regarded as 
a unique special case (linear dynamical system) to which our imaginary time trans- 
formation arguments may be applied. Fortunately, we can give a number of nonlinear 
models whose (Euclidean) inversion can be consistently implemented, [21]. The Eu- 
clidean connection goes beyond the confining vs scattering potential idea of ours and 
extends to periodic potentials as well. Examples from the instanton physics: static 
localized (kink) solutions of the 4 nonlinear field theory in one space dimension may 
be interpreted as Euclidean time solutions of the double well potential problem; the 
sine- Gordon kink may be interpreted as a Euclidean time solution of a plane pendulum 
problem. 

5.3 Time duality via analytic continuation in time 

The above procedure gives clear hints on how to connect the dual classical wave theory 
evolutions, associated with the previously discussed Hamiltonians H~[. We recall that 
in addition to the continuity equation, we infer the dual Hamilton- Jacobi equations 
d t s + (l/2m)(Vs) 2 ± V — and that there holds d t p = — V • (pv) with v(x,t) = 
(l/m)Vs(x,t). 

In the adopted notational convention, we define the initial data sq(x) = — Sq(x) and 
introduce an "imaginary time" transformation : 

ip(x,t) = p 1//2 exp(is/2mD) — > ip(x,t) = ip_ iso (x, —it) = (71) 

P-is ( x i —tt) exp[is_j So (x, — it)/2mD] = ~p l l 2 (x,t) exp[— s(x, t)/2mD] . 

We note that lim t | is^ iS() (x, —it) = i(—is )(x,0) = s (x). 

An analogous procedure for an analytic continuation in time has been worked out 
in the general context of the Euclidean quantum mechanics in Ref . [13] , however with 
no mention of its extension to the classically inspired Hamilton- Jacobi equation. 

Let us denote v = (1/m) Vs. Accordingly, the transformations implemented by (17TT) 
replace = J dxp[(m/2)v 2 + V] by — H ^ = J dxp[—(m/2)v 2 + V], with the very 
same function V(x) in both expressions. 

Clearly: 

d t p = -V • (pv) — > d t p = +V • (pv) (72) 
which is an obvious indicative of the time reflected (backwards) evolution. Analogously 

d t s + (l/2m)(Vs) 2 ± V = — ► d t s - (l/2m)(Vs) 2 + V = (73) 



17 



where, the time reflection t — > —t induces an expected form of the dual Hamilton- Jacobi 
equation: 

d t s + (l/2m)(Vs) 2 - V = . (74) 

The above discussed analytic continuation in time directly extends to the general 
pair of dual Hamiltonians, c.f. Section IV of Ref. [13]. The description be- 
comes even more straightforward, because in this case we connect pairs of linear par- 
tial differential equations. If ip{x,t) actually is a solution of the Schr"odinger equation 
i(2mD)d t = Hip, then 

i(j-i S0 (x, —it) = p^ 2 (x, t) exp[— s(x, t)/2mD] = 9*(x,t) (75) 

solves a backwards diffusion-type equation 

- (2mD)dA = H9, (76) 

while 

9(x,t) = p 1/2 (x,t)exp[+s(x,t)/2mD] (77) 
solves the forward equation 

(2mD)d t 9 = HQ . (78) 

In the above one may obviously identify D = h/2m. 

The whole procedure can inverted and we can trace back a non-dissipative quantum 
dynamics pattern which stays in affinity (duality) with a given dissipative dynamics, 
c.f. O. 



5.4 Diffusion- type processes 
5.4.1 Smoluchowski process 

The Hamiltonian appropriate for the description of dissipative processes (strictly speak- 
ing, diffusion-type stochastic processes) has the form 

H~ = J dxp [(m/2)v 2 - V - [m/2)u 2 ] (79) 

with the a priori chosen, continuous and bounded from below potential V(x). It is the 
functional form of V(x) which determines local characteristics of the diffusion process, 

Once the Fokker-Planck equation is inferred 

d t p = DAp - V(6 • p) , (80) 

where po(x) stands for the initial condition, we adopt the forward drift b = f/rwy of 
the process in the standard Smoluchowski form, characteristic for the Brownian motion 
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in an external force field f(x) = — VV. Here, 7 is a friction (damping) parameter and, 
instead of D = H/2m, we prefer to think in terms of D — k B T , /m / ~f where T stands for 
an (equilibrium) temperature of the reservoir. 

An admissible form of V — > / = — VV must be compatible with the Riccatti-type 
equation, provided the potential function V(x) has been a priori chosen: 



V(x) = m 



1 / f \ 2 ( f 



2 V m 7 J V m 7 



(81) 



The Fokker-Planck equation can rewritten as a continuity equation d t p = —V • j 
with the diffusion current j in the form: 

j = pv = -P- [f - k B TV In p] = ^ Vs . (82) 
m7 m 

We recall the general definition of the current velocity v = (l/m)Vs. 

Since the time-independent s = s(x) is here admissible, we have actually determined 

s = — — (V + ksTlnp) (83) 

7 

whose negative mean value F — —(s) determines for the Helmholtz free energy of the 
random motion, as follows: 

\U = 1 F = U-TS, (84) 

where S = ks S stands for the Gibbs-Shannon entropy of the continuous probability 
distribution, while an internal energy reads U = (V). 

Since we assume p and pVv to vanish at the integration volume boundaries, we get 

* = _( m7 ) (v 2 ) = -k B T(S) int < . (85) 

Clearly, the Helmholtz free energy ^ decreases as a function of time, or remains con- 
stant. 

The Shannon entropy S(t) = —(In p) typically is not a conserved quantity. We im- 
pose boundary restrictions that p, vp, bp vanish at spatial infinities or other integration 
interval borders and consider: 



which rewrites as follows 



where 



DS = (v 2 ) - (b-v) . (86) 

S = (S) int + (S) ext (87) 
k B T(S) mt = m 7 (v 2 ) > (88) 
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stands for the entropy production rate, while 

k B T(S) ext = - / / • jdx = -my (b ■ v) (89) 



(as long as negative which is not a must) may be interpreted as the heat dissipation 
rate:— J f ■ j dx. 

Let us consider the stationary regime S = associated with an (a priori assumed 
to exist) invariant density p*. Then, 



and 



Hence 



u = D V In p* 



- (l/k B T) VV = V In p* =>- p* = i exp[-V/£; B T] . (90) 



- 7S* = V + fc B T In p* =^ = - k B T \nZ = 7 F, (91) 

with Z = J exp(— V /ksT)dx. ^ stands for a minimum of the time-dependent Helmholtz 
free energy Because of 

Z = cxp(-^*/k B T) (92) 

we have 

p* = exp[(*„ - l/)/A; B T] . (93) 

Therefore, the conditional Kullback-Leibler entropy 7i c , of the density p relative 
to an equilibrium (stationary) density p* acquires the form 

k B TH c = -k B T / pln(— )dx = (94) 

In view of the concavity property of the function f(w) = —wlnw, TC C takes only 
negative values, with a maximum at 0. We have ^* < ^ and k B TTi, c = — \I/ > 0. TC c (t) 
is bound to grow monotonically towards 0, while *&(t) drops down to its minimum ty* 
which is reached upon p*. 

The Helmholtz free energy minimum remains divorced from any extremal prop- 
erty of the Gibbs-Shannon entropy. Only the Kullback-Leibler entropy shows up an 
expected (growth) asymptotic behavior. See e.g. also [2]. 

Note that properties of the free Brownian motion can be easily inferred by setting 
b = in the above discussion. Then, the diffusive dynamics is sweeping and there 
is no asymptotic invariant density, nor a finite minimum for \l/(t) which decreases 
indefinitely. 
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5.4.2 Reintroducing duality 

To set a connection with the previous time duality ("imaginary time" transformation) 
framework, we need only to observe some classic properties of Smoluchowski diffusion 
processes. Once we set b = — 2Z)V$ with $ = a substitution: 

p(x,t) =^(z,t)exp[-$(x)] (95) 

with 9* and $ being real functions, converts the Fokker-Planck equation into a gener- 
alized diffusion equation for 6**: 

dtO* = DA9* - (96) 
ImD 

and its (here trivialized in view of the time-independence of $) time adjoint 

d t 9 = -DA6 + ^-9 (97) 
for a real function 8(x,t) = exp[— $(x)], where 

si = ^ + v ' t ) = i, i( v *) 2 - A *]- (98) 

Let us note an obvious factorization property for the Fokker-Planck probability density: 

p(x,t) = 9(x,t) ■ 0*(x,t) (99) 

which stays in affinity with a quantum mechanical factorization formula p = 
albeit presently realized in terms of two real functions 9 and 9*, instead of a complex 
conjugate pair. In view of (at this point we restore an original notation of Section 4): 

p 1/2 {x,t) exp[-s{x,t)/2mD] = 9*(x,t) (100) 

we immediately recover 

s = (2mD)[<$> - (1/2) hip] (101) 

in conformity with the previous definition Eq. (83). If there are no external forces, 
$ disappears and we are left with the free Brownian motion associated with s = 
—mD In ~p. 

For the record it is useful to mention explicit transformations between Green's func- 
tions appropriate for quantum motion and transition probability densities of standard 
diffusion type processes, [25]. Explicit examples of the free dual dynamics and those 
in the harmonic potential, have been worked out there. 
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6 Conclusions 



An analytic continuation in time (or an "imaginary time" transformation) stands for 
a mapping between two different types of dynamics, both running with respect to the 
equally " real" (istic) time. In our case, the modular Schrodinger equation has been 
a unifying departure point for an analysis of, sometimes not quite expected, affinities 
between dynamical patterns of behavior generated by the same primary nonlinear field, 
but typically considered disjointly - as research problems on their own, for properly 
selected coupling parameter values. 

One may possibly take the view that k = 0, 1, 2 cases are just formally analo- 
gous mathematical descriptions of different physical systems. Our standpoint is that 
a primary dynamical system is a modular Schrodinger equation with an arbitrarily 
adjustable coupling parameter. Different coupling regimes refer to physically different 
patterns of behavior, but there is a deep intertwine between them, to be further ex- 
plored. The global changes of properties of solutions of nonlinear dynamical equations 
as the control parameter is varied are a routine wisdom for nonlinear system experts. 
This property is shared by the modular Schrodinger equation as well. 

If one takes seriously the dynamical duality (or time duality) concept, the models 
considered here should not be viewed independently anymore. One can trace the 
dynamical patterns of one model in terms of those for another, and in reverse. Even, 
if at the first glance, the pertinent dynamics patterns may seem to have nothing in 
common. 

It is the Hamiltonian analysis of the nonlinear field (e.q. Schrodinger equation) 
which has led to a variety of emergent Hamiltonian motion scenarios. Two classes of 
them, we call them confining or scattering, stay in close affinity on quantum, classical 
and stochastic (dissipative) dynamics levels of description. 

Euclidean methods are often used in various (especially quantum) branches of the 
statistical physics research of equilibrium and near-equilibrium phenomena. Our dis- 
cussion was basically concentrated on the real time flow notion which, even after a 
Euclidean (imaginary time) transformation, still remains a real(istic) time flow, albeit 
with a new physical meaning. 

It is worth mentioning that an an independent, quantum theory motivated approach 
[26] (deformation quantization with an "imaginary transformed" deformation constant 
h) has obvious links with the Euclidean map viewpoint towards diffusion-type processes, 
explored in the present paper. 

We would like to point out that, quite aside of existing and prospective physical 
implementations, the sign-inversion issue for the conservative potentials has a deeper 
mathematical meaning whose role we have slightly diminished, not to overburden the 
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text with a strong admixture of an advanced functional analysis. The dynamical semi- 
group indications in the present paper were intended to tell the mathematically in- 
terested reader under what circumstances one can be sure of the existence (modulo 
suitable time interval limitations) of solutions that are connected by an imaginary 
time transformation. The semigroup notion cannot be hastily extended to the classical 
dynamics, nonetheless an imaginary time link still works there. 

Acknowledgement: I am indebted to Leon Brenig for pointing out to my atten- 
tion his paper [1] and enlightening correspondence. 
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